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We discuss the dynamics of a Bose-Einstein condensate 
during its nondestructive imaging. A generalized Lindblad 
superoperator in the condensate master equation is used to 
include the effect of the measurement. A continuous imag- 
ing with a sufficiently high laser intensity progressively drives 
the quantum state of the condensate into number squeezed 
states. Observable consequences of such a measurement- 
induced squeezing are discussed. 
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Since its birth, quantum mechanics has led to an inter- 
pretational debate on the role played by the measurement 
process in its structure and its relationship to classical 
mechanics developed for macroscopic systems This 
debate has been enriched by the realization of new exper- 
imental techniques spanning from quantum jumps in sin- 
gle ion traps to macroscopic entangled states in various 
quantum systems. Recently, the production of atomic 
Bose-Einstein condensates of dilute atomic gases has also 
paved the way to the study of dynamical phenomena of 
macroscopic quantum systems with the precision charac- 
teristic of atomic physics || . 

Two optical techniques, absorption and dispersive 
imagings, have been used to monitor the dynamics of 
a Bose-Einstein condensate 0]. In the former the con- 
densate interacts with a light beam resonant (or close to 
resonance) with an atomic transition. The output beam 
is attenuated proportionally to the column density of the 
condensate - the condensate density integrated along the 
line of sight of the imaging beam. The absorption of 
photons heats the condensate then strongly perturbing 
it, and in general a new replica of the condensate has to 
be produced to further study its dynamics. This mea- 
surement is, in the language of quantum measurement 
theory, of type-II since it destroys the state of the ob- 
served system and forbids the study of the dynamics of 
a single quantum system || . Repeated measurements on 
a Bose-Einstein condensate or, at the limit, its continu- 
ous monitoring are instead possible using its dispersive 
features, for instance through phase-contrast [Q or in- 
terference H imaging techniques. Off-resonance light is 
scattered by the condensate which induces phase-shifts 
thereby converted into light intensity modulations by ho- 
modyne or heterodyne detection. The off-resonant na- 



ture of the atom-photon interaction allows for a very low 
absorption rate and therefore low heating of the conden- 
sate. Thus, multiple shots of the same condensate can 
be taken - a type-I measurement - allowing to study with 
high accuracy several phenomena, like its formation in 
non-adiabatic conditions ||, short and long wavelength 
collective excitations Q, vortices and superfiuid dynam- 
ics [gj . The effect of the measurement process is typically 
neglected in these analysis. A first attempt to include the 
measurement process in a two-mode configuration for the 
condensate has been discussed in ||. Our main goal is 
to include the atom-photon interaction process present in 
dispersive imaging into the intrinsic dynamics of the con- 
densate. We show that the measurement process induces 
number squeezing of the quantum state of the conden- 
sate, and we provide realistic estimates that suggest that 
the phenomenon could be observed with current experi- 
mental techniques. 

Let us start the analysis with the effective interaction 
Hamiltonian between the off-resonant photons and the 
atoms, written as: 
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where h is the condensate density operator, and E the 
electric field due to the intensity I of the incoming light. 
The coefficient xo represents the effective electric sus- 
ceptibility of the atoms defined as xo = X 3 S/2tt 2 (1 + 5 2 ), 
where we have introduced the light wavelength A and 
the light detuning measured in half-linewidths T/2 of the 
atomic transition, S = (w — w at )/(T/2). Equation (|I|) 
allows us to write the reduced master equation for the 
atomic degrees of freedom by a standard technique, i.e. 
by tracing out the photon degrees of freedom p0[ . The 
photons are assumed to be in a plane-wave state with 
momentum along the impinging direction, correspond- 
ing to a wavevector k = koz orthogonal to the imaging 
plane x — y. Unless tomographic techniques are used, the 
image results from a projection of the condensate onto 
the x — y plane, by integrating along the z direction. 
This demands to project the dynamics of the conden- 
sate into the imaging plane. In order to write a closed 
2-D master equation to describe the x — y dynamics we 
assume the condensate wavefunction to be factorizable 
as ip(x,y,z) — 4>(x , y) A(z) . Such factorization holds 
if the confinement in the ^-direction is strong enough 
to make the corresponding mean-field energy negligible 
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with respect to the energy quanta of the confinement, 
i.e. hto z » Ng, as recently demonstrated in O], where 
g = AirTi 2 a/m, with a the s-wave scattering length, and 
to z the angular frequency of the confinement harmonic 
potential along the z direction. The resulting reduced 
master equation in the imaging plane is written as p^ ]: 



d 2 n f d 2 r 2 K(r 1 - r 2 )[n(n), [n(r 2 ), p}}, (2) 



dp 
~dt 



where ft(r) = ^^(r) is the 2-D density operator. The 
effect of the measurement is taken into account through 
the second term in the right-hand side of Eq.(|^). This 
equation preserves the total number of atoms, and cor- 
responds to a quantum nondemolition p3 coupling be- 
tween the atom and the optical fields p.|f 4| 1G| . The mea- 
surement kernel K has the expression: 

K(r) = J d 2 kcM-ek 4 /4k 2 ) exp(zkr), (3) 

where £ is the lengthscale of the condensate in the z 
direction, the width of the Gaussian state A(z) under 
the abovementioned approximation. Equation (^|) holds 
for a condensate having thickness in the z direction 
£ ^> A. If the measurement kernel were a local one, 
K{t\ — r 2 ) ~ <5(ri — r 2 ), Eq.(||) would reduce to a Lind- 
blad equation for the measurement of an infinite num- 
ber of densities n(r). This assumes that no spatial cor- 
relation is estabilished by the photon detection. How- 
ever, the ultimate resolution limit in the imaging sys- 
tem depends on the photon wavelength, regardless of the 
pixel density of the detecting camera. The resolution 
lengthscale follows from Eq.(H) as a width of the kernel 
Ar = (27r 2 ^//co) 1 ^ 2 = (^C^) 1 j t ne geometrical average 
of the light wavelength and the condensate thickness £. 

When the measurement outweighs the self dynamics, 
the wavefunction of the condensate is driven towards an 
eigenstate of the measurement apparatus. In that case 
the master equation (||) is solved by the eigenstates of 
the density operator h(r), which are then pointer states 
0, * f (ri) • ■ • ^(rjv)^). This state describes a "gas" of 
N atoms localized at the points ri . . . rj\r. These states 
are manifestly not mean-field states so they cannot be 
described by any generalization of the Gross-Pitaevskii 
equation. A measurement strong in comparison to the 
Hamiltonian H2B will project an initial mean field con- 
densate state onto one of the pointer states with a prob- 
ability distribution for different r's given by the initial 
mean field wave function p(r) = (j)*<f>{y). However, this 
dramatic change of the condensate state will not directly 
affect the condensate image since the homodyne current 
is given by the convolution /(r) ~ J d 2 r'K(r — r')(ft(r')), 
which coarse-grains the density distribution over the res- 
olution of the kernel Ar 2 . The only effect on the image 



will be small statistical fluctuations: if the average num- 
ber of atoms within a given resolution area of Ar 2 is n, 
then after the localization the number will be in the range 
n ± y/n. 

When the system Hamiltonian H2B is strong enough, 
it can change pointer states |jl8). To quantify this phe- 
nomenon and the competition between localization in- 
duced by the measurement and the derealization due to 
i? 2 D we discretize Eq. (|2j) introducing a 2-D lattice with 
the natural choice for the lattice constant set by the res- 
olution lengthscale Ar, 



dp 



= " [ _ n ^^ih + V,p}-Sj2[ni,[ni,p\]. (4) 



(k,l) 



Here ^1 is an annihilation operator and hi = ^\^i 
is the number operator at a lattice site The fre- 
quency of hopping between any nearest neighbor sites 
(k,l) is u) as h/2mAr 2 , which is TT 1 times the charac- 
teristic kinetic energy. The potential energy operator is 
V = Y^i Uihi + Ghf, where Ui is the trapping potential 
and G — g/£Ar 2 . The effective measurement strength is 
S « Jd 2 rK(r)/Ar 2 = (2tt / ' Ar) 2 {n X lk I /2hc). 

The Lindblad term in Eq. (^) is a sum of Lindblad 
terms for different lattice sites. Its pointer states are 
the Fock states \{Ni}) = \Ni,N 2 ,...) with definite oc- 
cupation numbers Ni at every site I. An initial pure 
mean- field state can be expanded in the Fock basis, 
|-0(O)> = E{jv,} i>{Ni}(0)\{Nl})- If thc initial mean-field 
wavefunction of the condensate is , then the initial ex- 
pectation value of the number of atoms at a site I is 
ni(0) = N<p*(pi. Each site has a finite dispersion in the 
number of atoms, <t;(0) ~ [rii(O)] 1 / 2 . For large ni(0) the 
distribution of Ni can be approximated by a Gaussian. If 
there were no self-Hamiltonian the Lindblad term would 
drive the state of the system into a mixed state of dif- 
ferent Fock states |{JVj}). The mixed state describes an 
ensemble of outcomes of different realizations of the ex- 
periment. In a given realisation the outcome would be a 
definite Fock state |{ATj}) with random iV; in the range 
ni(0)±<Ji(0). To describe a single realisation of the exper- 
iment we use a (Stratonovich) unraveling of the master 
equation by a stochastic Schrodinger equation (SSE) JnJ 



-Tp{N l} Y,[- S ( N '- n rf + S(T ? + ( N l- n ^] > ( 5 ) 







where the homodyne noises have averages 0i(t) 
and e^Jj^hJti) = 2SS hM S(t 1 - t 2 ). Here n x = 
W{jv,}| 2 , <Ji = T, {Nl} (Ni - nif\i, {Nl] \ 2 , h is 
the matrix element of the hopping Hamiltonian and 
V {Nl} = J2ii u i N i + GN i) is the potential energy. The 
unravelling does not change the pointer states [|20(|. 
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for a 1-D 3-mode configuration with 
S = 2s- 1 , ni(0) = 100 and of(0) = 66.6. We show four 
single realizations of the experiment for different hopping 
frequencies. The thick solid line is the analytic solution 
of (*) = ff? (0)/[l + 45^(0)] for w = s" 1 . 

Equation (||) can be solved provided that we neglect 
the hopping term, h = 0. In this case Fock states are 
fixed points because they are characterized by Ni = ni 
and a i = 0. SSE is satisfied exactly by the ansatz 



^{iVi}(*) 



n 



[2W(*)l 1/4! 



(6) 



when the parameters satisfy the equations of motion 



dt 



(7) 



and <P{Ni} — J2i NiArg(4>i) are constant phases. The 
ansatz remains normalized as long as of 3> 1 when sum- 
mations over TV; can be replaced by integrals. For an 
initial of (0) ~ «;(0) the dispersions shrink to zero (a 
Fock state) as af{t) = of (0)/[l + 4of (0)Si]. Intuitively, 
the continuous reduction of number fluctuations is due 
to the weak nature of the measurement, and takes place 
on a time scale ~ 1/Sof (0) 



While of 's shrink down to 



of ~ 1 the means ni make random walks in the range 
ni(0) ± 07(0). When of -C 1 our ansatz breaks down but 
the numerical simulations for 2- and 3-site models show 
that ni becomes frozen. 

For non zero h, localization in a Fock state may be 
inhibited due to the hopping between nearest neighbor 
sites. Here a scaling argument is used to estimate the 
dependence on the physical parameters of of (oo), the dis- 
persion at which the measurement term and the hopping 
term in the SSE balance each other. We assume that the 
density distribution is smooth, i.e. any two nearest neigh- 
bor sites k,l have close occupation numbers, riy. ps raj. 
For large occupation numbers, the element of the hop- 
ping Hamiltonian between k and I can be approximated 



by ujy/NkNi m Lo^/n^ni « wrij. The hopping mixes the 
amplitudes ip{N t } an d ip{N'} which have negligibly close 
potential energies, Vjjv,} ~ V{N[}- We again assume that 
of ^ 1 so that we can treat the occupation numbers Ni 
as continuous variables. After neglecting the noise term, 
the SSE can be written as an eigenvalue (Q) problem 



ifiV>{jv,} = ^{w,} ^2 [~S(Ni - ni) 2 - 
i 

~d 2 ip {Nl} d 2 tfj {Nl} 



Saf 



-iu y/n k ni 

(k,l) 



d 2 ^ { 



Ni} 



dN 2 



dN k 8N t 



(8) 



This equation is made parameter independent by rescal- 
ing Ni = Niiurn/S) 1 / 4 and ft = ft^SuTOj) 1 / 2 , so the 
stationary solution of the rescaled equation must scale 
as 07(00) ~ (ujni/ S) 1 / 4 . In Figures 1 and 2 we show 
results of exact numerical simulations on a 1-D 3-site pe- 
riodic lattice with 300 atoms which confirm the validity 
of these predictions. The solution without the hopping 
term, af(t) = of(0)/[l + 4Sfof(0)], reaches its asymp- 
totic value 07(00) after a time 77 — [Su;raz(0)] _1 / 2 . 

We estimate the effect assuming the following param- 
eters for the condensate and its imaging, relevant for the 
case of 87 Rb: m = 1.4 10~ 25 kg, a = 5.8nm, A = 780nm, 
Xo = 10 _23 m 3 , laser intensity I = 100mW/cm 2 , and a 
total number of atoms N = 10 7 . The width of the Gaus- 
sian A(z) is assumed to be £ = 10/^m. The resolution of 
the kernel (||) is Ar = (irXt;) 1 / 2 = 5um. We also get uj — 
15s" 1 , S = 65s" 1 , and G/Ti = 0.2s- 1 < S. We assume 
that the initial condensate has a size of (50/xm) 2 which 
implies 100 occupied lattice sites with nj(0) » 10 7 /100 = 
10 5 atoms per site and 07(0) ~ [n^O)] 1 / 2 ~ 300 atoms per 
site. Due to the measurement the dispersion shrinks 25 
times down to 07 it = 00) = (wn;(0)/ S) 1 ^ 4 w 12 atoms 
per site in a timescale 17 = (S'wn;(0)) _1 / 2 « 100/is. Until 
now we have neglected the off-resonance Rayleigh scat- 
tering due to the measurement. The same methods lead- 
ing to Eqs.(|],[|) also predict a Rayleigh depletion time of 
tjj ~ 7ms, i. e. 70 times longer than 77 . 

The dramatic 25-fold shrinking of the dispersion in oc- 
cupation number results in a quasi-Fock state. An exact 
Fock state would be . . . &j |0), where each atom can 
be assigned to exactly one lattice site and as such it is 
localized within the kernel resolution of Ar 2 . The local- 
ization in an exact Fock state would cost a lot of kinetic 
energy. In our guasi-Fock state, as shown by the calcula- 
tion of the expectation value of the hopping term in our 
ansatz state (||), as long as all uj > 1 the energy remains 
very close to the energy of the initial mean-field state. 

A direct experimental test of our predictions and in 
particular of a simplified lattice model (^J) can be made 
in the 1-D array of weakly linked mesoscopic traps of 
the experiment described in Ref. pl[ . In this experi- 
ment the traps are created by a standing optical wave 
of wavelength A s t = 840nm superimposed on a confin- 
ing harmonic potential due to a magnetic trap. Given 
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a transverse trapping angular frequency of 2ir x 120s 1 
we estimate the transverse length of the condensates in 
each trap as the corresponding oscillator length, equal 
to £ = lfim. In the experiment each condensate has a 
1/e width of Al — A s t/6, therefore we estimate G/Ti — 
g / (2TiAl(2t;) 2 ) w 40s -1 . Assuming one is in the range 
of parameters of the experiment for which the ground 
state of the condensate is close to a mean-field (quasi- 
coherent) state (for instance squeezing s ~ 3 in Fig- 
ure 2 of Ref. ^J), and that in each trap there are 
n;(0) « 100 atoms of 87 Rb, the hopping frequency is 
lu = Gni(0)/hs 2 w 500s- 1 . 



condensate. We have shown that dispersive imaging 
with a sufficiently high laser intensity results in number 
squeezing of the condensate state. Our prediction can 
be tested in present-day experiments and could result 
in a significant improvement toward the implementation 
of Heisenberg-limited atom interferometry 0] , and in a 
more general understanding of the back-action of other 
quantum- limited devices [ p^[ . 
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FIG. 2. Asymptotic variance a 2 (oo) versus (uini(0) / S) 1 ' 2 
for a 1-D 3-mode configuration, with S = 2s _1 and 
n((0) = 100. The scaling of(oo) ~ lj 1//2 applies for 
(ujm(0)/S) 1/2 < erf(0) = 67 when u < Saf(0)/ni. The 
straight line is a guide for the eye. In the opposite case of (oo) 
does not increase above the dispersion of the initial mean-field 
state. The of (oo) are averages over the noisy values at large 
t of trajectories in Fig.l. 

Suppose now that we add phase contrast imaging to 
the above experimental setup. The kernel resolution is 
Ar = {tt^X) 1 / 2 « 1.5^tm, which implies there are 2 op- 
tical traps within the kernel resolution. The effective 
number of atoms within Ar is therefore n pc i = 200, the 
effective self coupling is G pc i/h = G/2h = 20s -1 and 
the effective hopping frequency is w pc i = lo/2 = 250s -1 . 
For a laser intensity of, say, / = 100mW/cm 2 , we obtain 
S = 500s _1 and c;(oo) = 3 atoms per site after a local- 
ization time of r; = 200/zs but much before a depletion 
time of td = 7ms. The final dispersion in the number 
of atoms per optical trap is cr(oo) = (aiipo) 2 /2) 1 / 2 = 2. 
The final squeezing factor will be ct 2 (0)/ct 2 (oo) ss 50, 
where we take of(0) = n/(0). Such an enhancement of 
squeezing due to the measurement should be observable 
with the techniques of the experiment [ pl| 

In conclusion, we have applied the theory of open quan- 
tum systems to include the back-action of a nondemoli- 
tive measurement into the dynamics of a Bose-Einstein 
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